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The detection of non-baryonic dark matter through its gamma-ray annihilation in the centre of our galaxy 
has been studied. The gamma fluxes according to different models have been simulated and compared to those 
expected to be observed with the Alpha Magnetic Spectrometer (AMS), during a long-term mission on board of 
the International Space Station. Under the assumption that the dark matter is composed of the lightest, stable 
supersymmetric particle, the neutralino, the results of the simulations in the framework of mSUGRA models, 
show that with a cuspy dark matter halo proflle or a clumpy halo, the annihilation gamma-ray signal would 
be detected by AMS. More optimistic perspectives are obtained with the Anomaly Mediated Supersymmetry 
Breaking (AMSB) model. The latter leads also to a cosmologically important ®Li abundance. Finally, the 
discovery potential for the massive Kaluza-Klein dark matter candidates has been evaluated and their detection 
looks feasible. 



1. Introduction 

The nature of dark matter is one of the out- 
standing questions and challenges in cosmology. 
The existence of cosmological dark matter is re- 
quired by a multitude of observations and ar- 
guments, such as the excessive peculiar veloc- 
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ities of galaxies within clusters of galaxies, or 
the observations of gravitational arcs indicating 
much deeper gravitational potentials within clus- 
ters than those inferred by the presence of the 
luminous matter jll2j . On the galactic scale, ex- 
tensive dark matter halos are required to explain 
the observed rotation curves in spiral galaxies, or 
the velocity dispersion in elliptical galaxies |3I4| . 
Furthermore, Big Bang nucleosynthesis predicts 
a fractional contribution of baryons to the criti- 
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cal density, fif,, significantly smaller than the total 
fljn in form of clumpy matter. The Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) has provided 
the most detailed measurements of the Cosmic 
Microwave Background (CMB) anisotropics [5]. 
In the framework of the Standard Cosmological 
Model, WMAP quotes a total matter density of 
rim= 0.27±0.04 and a baryon density of fib — 
0.044±0.004, which confirms that most of the 
matter is non-baryonic, in agreement with the 
results obtained from primordial nucleosynthesis 
studies. 

Various non-baryonic dark matter candidates 
require physics beyond the standard model of par- 
ticle physics (for a recent review see e.g. N- 
body simulations of structure formation [7] sug- 
gest a non-relativistic. Weakly Interacting Mas- 
sive Particle (WIMP) as a dark matter compo- 
nent, thus favoring the Cold Dark Matter sce- 
nario IHEl- The WMAP measurement of the 
density of the non-baryonic dark matter provides 
constraints in the range of 0.095 < ftcDMh^ < 
0.129, at the 2 a level. 

Supersymmetric theories offer an excellent 
WIMP candidate, which satisfies the CDM 
paradigm and the constraints on flcDivh namely 
the neutralino (xi) of the Minimal Supersymmet- 
ric Standard Model (MSSM), assumed to be the 
Lightest Supersymmetric Particle (LSP) and sta- 
ble due to i?-parity conservation ^HI- At present, 
lower limits on the LSP neutralino mass in the 
MSSM are about 50 GeV from LEP experiments 
(although the mass may be significantly smaller 
depending on the assumptions relative to gaug- 
ino mass universality). Less conventional scenar- 
ios than the neutralino within the minimal super- 
gravity (mSUGRA) context have been proposed 

• In the Anomaly Mediated Supersymmetry 
Breaking (AMSB) scenario il2|il3j the neu- 
tralino LSP is predominantly a Wino [the 
supersymmetric partner of the electrically 
neutral component of the SU{2)l gauge 
bosons]. Endowed with a relatively large 
annihilation cross section this particle may 
constitute the bulk of the dark matter when 
subsequently to its thermal freeze-out, it is 



further generated nonthermally (e.g. via Q- 
ball evaporation or gravitino decay). By 
virtue of its large annihilation cross section 
the wino may lead to possibly large gamma 
fluxes HH. 

• In extra-dimension models, ultimately mo- 
tivated by string theories, it has been ar- 
gued ^S] that the lightest Kaluza-Klein ex- 
citation can provide under certain condi- 
tions, a very good CDM candidate. In the 
present paper we will restrict ourselves to 
the possibility of low scale extra-dimensions 
as an extension of the non-supersymmetric 
standard model IB' with a perfectly viable 
dark matter Kaluza-Klein particle jT2] . 

In this paper, we present the predicted 7- 
ray fluxes from the Galactic Centre from neu- 
tralino annihilations in the frame of mSU GRA 
and AMSB models, as well as from Kaluza-Klein 
dark matter annihilations. The predicted fluxes 
are used to assess, for the different scenarios, the 
discovery potential for non-baryonic dark matter 
provided by a three-year observation of the dif- 
fuse 7-ray differential spectrum by the AMS on 
the ISS. 

2. Model descriptions and simulations 

The limited knowledge of dark matter structure 
and the density profile near the Galactic Center 
represent the principal astrophysical uncertain- 
ties when evaluating the discovery potential of 
the dark matter through indirect detection. On 
the other hand, the predicted 7-ray fluxes depend 
on the assumptions made within the framework 
of the particle physics models associated with the 
different dark matter candidates. 

2.1. Dark matter halo parametrization 

One may parameterize the mass density profile 
of our Galaxy by the following equation: 

assuming a simple spherical Galactic halo. An 
isothermal profile with core radius a corresponds 
to 7 = 0, a = 2 and e = 1 as proposed by 
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[TH| . A Navarro, Frenk and White (NFW) profile 
[TU| is obtained with 7 = 1, a = 1, and e = 2, 
whereas Moore's distribution 20 is recovered if 
7 = e = 3/2 and a = 1. Only the NFW and 
Moore models are considered in this study. The 
two models predict large values of the neutralino 
density in the Galactic Centre (GC). 
The parameters of the halo modeling are: 

• Rq - distance from Earth to GC, 

• Po - halo density at Rq, 

• a - the core radius - for r < a the halo 
density is constant and equal to p{a) in case 
of the isothermal parametrization. 

As shown in |22| . po and a cannot be chosen 
arbitrarily. The total mass of the Galaxy re- 
stricts the {po, a) parameter space. For the NFW- 
standard model the generic parameter values are: 
Rq = 8.0 kpc, Pa = 0.3 GeV/cm^, a = 20 kpc. 
Another possible combinations of {po, a) param- 
eters allow, given the uncertainties: Rq — 8.5 kpc, 
Po = 0.4 GeV/cm^, a = 4 kpc (NFW-cuspy). The 
two configurations were considered. The values 
for the Moore profile have been chosen as follows: 
Rq = 8.0 kpc, Po = 0.3 GeV/cm^, a = 28 kpc. 

The WIMPs located aroimd the Galactic center 
should annihilate and produce high energy pho- 
tons. The corresponding photon flux near the 
Earth, <i>^ - per unit time, surface, and solid an- 
gle - may be expressed as 

1 (av) N^, f ■? / N , /„x 

4^ ^ ^wimp J los 

In equation 2, m^imp is the mass of the WIMP- 
dark-matter candidate; (av) denotes the ther- 
mally averaged annihilation rate ; pwimp (r) is the 
mass density of the dark matter and r is the dis- 
tance from the Galactic center. The flux $^ is 
proportional to the number of annihilations per 
unit time and volume, {av) P^impi^) / ''^wimp ^^'^ 
to the number of secondary photons per annihila- 
tion, A''^. Finally to obtain the flux at the Earth 
it is necessary to integrate the WIMP density 
squared along the line-of-sight (los) connecting 



the observer to the Galactic Center. The inte- 
gral can be expressed in the form: 



J(i?) = 2 y p^ iJs^+R^j ds , (3) 

assuming a spherical halo with radial extension 
i?o- The coordinate s extends along the line of 
sight. R is the radial distance from the centre of 
the Galaxy for such a direction. 

Because the density decreases steeply at large 
distances (see Eq. ^ , a radial cut-off Rc has been 
applied. It has been checked that our results are 
not sensitive to the value of Rc, which is set to 8.0 
or 8.5 kpc depending on the chosen halo profile. 

We integrate the function J over a solid an- 
gle around the Galactic center, subtended by the 
detector acceptance, e.g. a circular region with 
angular radius ^obs: 

/■^obs 

E = 27r / J{R) ain0 d0 , (4) 
Jo 

where R/Ro — tanf? ~ 9. Thus the resulting 
value for Ij - flux of high-energy photons col- 
lected per unit of time and surface - can be writ- 
ten as: 

— ^3.98 X 10 photons cm s ^) X 



where S19 denotes S expressed in units of 
10^^ GeV^ cm^*^ and integrated over the AMS ac- 
ceptance. Ail — 10~^ sr. We have integrated the 
relation 10} as a function of galactic halo profiles 
leading to the following results: 

• For a NFW-standard profile : S19 = 2.7 
10^, thus J{0){An) = 1.2 10^ 

• For a NFW-cuspy profile: E19 = 117.7 10^ 
thus J(0)(Arj) = 50.0 10^ 
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• For a Moore profile: Eig = 336.7 10^, thus 
J(0)(AO) = 142.9 10^ 

where J{Q){AU) corresponds to the notation used 
in |22| . For integration, we have defined an inner 
cutoff radius at Rc = 10~^ pc such that for R > 
Rc the flux vanishes. 

For completeness, two astrophysical factors 
may enhance the expected gamma fluxes from 
neutralino annihilations: 

• dumpiness of the dark matter halo as indi- 
cated by N-body simulations 

• the presence of a Supermassive Black Hole 
(SBH) with a mass of ~ 2.6 x 10^ Mq creat- 
ing unstable conditions due to baryon infall 
by the adiabatic compression process, stud- 
ied by ^ 

The overall enhancement factor of the expected 
flux is estimated in |26 | to be between 5 and 
100, depending on the dumpiness of the galactic 
halo. The enhancement of the annihilation signal 
in presence of a spike in the dark matter halo is 
significant with respect to ordinary dark matter 
cusp, even in case of gravitationnal scattering of 
stars and the self-annihilating dark matter parti- 
cles, as pointed by |23I24I25| . 

2.2. Models for WIMP candidates 

In this section, we describe briefly the physics 
models associated with selected dark matter can- 
didates, including the methodology and assump- 
tions used to investigate the different hypotheses. 

2.2.1. mSUGRA parameterization 

The two supersymmetric scenarios considered 
belong to the class of models where supersymme- 
try (SUSY) breaking is effectively communicated 
to the visible sector via (super)-gravitational ef- 
fects. We use the conventional mSUGRA sce- 
nario |11) with common values for the soft su- 
persymmetry breaking scalar and gaugino masses 
and trilinear couplings, mo,TOi/2,^Oj taken as 
initial conditions at a given high energy universal- 
ity scale. We require the three gauge couplings to 
take a common value at a uniflcation scale Mqjjt 
and, for simplicity, identify this scale with the 



universality scale of the soft SUSY parameters. 
With these initial conditions and with the value 
of tan /? (the ratio of the two Higgs vacuum expec- 
tation values ^jfi-^) deflned at the electroweak 
scale, the relevant low energy quantities are ob- 
tained throug the running of the parameters from 
Mqut down to a scale of the order of the elec- 
troweak scale. Electroweak symmetry breaking 
is then required at that scale with the correct Z 
boson physical mass, thus fixing the supersym- 
metric /i parameter (up to a sign) and its soft 
supersymmetry breaking counterpart. 

The flux predictions are obtained by use of 
computational tools which allow to scan over var- 
ious SUSY parameters. This is achieved through 
an interface of the two codes, DarkSUSY 122] 
and SUSPECT"^, which we dub hereafter DSS 
(DarkSUSY-SUSPECT). Signiflcant features of 
particle physics and cosmology are thus combined 
in our approach, taking into account various phe- 
nomenological constraints (consistency of the top, 
bottom and r masses, present experimental lim- 
its on the superpartner and Higgs masses, lim- 
its from h — > S7, no charged LSP, relic den- 
sity constraints) , some of which are implemented 
in SUSPECT and others in DarkSUSY. We have 
checked in the mSUGRA framework the compat- 
ibility of the results obtained with the DSS soft- 
ware package and the ISASUGRA interface pro- 
vided with DarkSUSY. 

2.2.2. Anomaly Mediated Supersymmetry 
Breaking parameterization 

The Anomaly Mediated Supersymmetry 
Breaking (AMSB) is a gravity- mediated mech- 
anism where the SUSY breaking is communi- 
cated to the observable sector by the super- Weyl 
anomaly |12I13| . In particular, the masses of 
gauginos are generated at the one-loop level as in 

Mr=h[^y<M> (6) 

where ai are the gauge coupling constants and bt 
the associated /3-function coefficients. M is the 
auxiliary field in the supergravity multiplet whose 
vacuum expectation value < M > is expected to 
be of the order of the gravitino mass 777,3/2, the 
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latter being generically in the range: 10 TeV < 
m3/2 < 100 TeV. 

In the minimal AMSB model, the sleptons suf- 
fer typically from a tachyonic problem. One way 
to fix this problem is to add a scalar mass param- 
eter mg, accounting for a non- anomalous contri- 
bution to the soft SUSY breaking. 

The most important message from the gaugino 
mass formula above is the hierarchy: 

Ml : M2 : M3 = 2.8 : 1 : 8.3 

as opposed to Af 1 : M2 : Mi = 1:2:7 which is 
expected for the gravity- or gauge-mediated mod- 
els. This implies that the lightest neutralino ( Xi) 
and the lightest chargino ( xf) ^re almost pure 
Winos and consequently mass-degenerate. 

2.2.3. Kaluza-Klein Dark Matter 

Models with compact extra dimensions predict 
several new states, the Kaluza-Klein (KK) exci- 
tations. In the case of Universal Extra Dimen- 
sions (UED) jT^, all Standard Model fields can 
propagate in the bulk, and their effective four- 
dimensional interactions with the KK states con- 
serve a quantum number associated with the lat- 
ter. The conservation of this KK number implies 
that the KK modes cannot decay exclusively into 
Standard Model particles; the lightest KK mode 
(LKP) wiU thus be stable |16I3H . Moreover, the 
LKP, when electrically neutral and with no bary- 
onic charge, provides a viable dark matter can- 
didate PZj. The mass of the LKP dark matter 
particle, like all other states of the KK tower, is 
inversely proportional to the compactification ra- 
dius R. Accelerator electroweak measurements 
constrain rather weakly the UED scenario, since 
in this scenario the observables are sensitive only 
to the virtual effects of the KK modes. The lower 
mass bound leads to i?-^ > 280 GeV I16j. The 
most promising LKP dark matter candidate is as- 
sociated with the first level of KK modes of the 
Hypercharge gauge boson B'-^^ . In our calculation 
we consider the relic density of B^^^ of ^\ lead- 
ing to a B^^^ lower bound mass bound constraint 
of 400 GeV. 



3. AMS gamma detection and sensitivity 
3.1. The AMS-02 experiment 




Figure 1. Schematic view of the AMS-02 ex- 
periment which will operate on the International 
Space Station. 



The main elements of AMS-02 detector |221 
which are shown in Figure ^ include: a super- 
conducting magnet, a gaseous transition radia- 
tion detector (TRD), a silicon tracker (Tracker), 
time-of-flight hodoscopes (TOF), a ring imag- 
ing Cerenkov detector (RICH), an electromag- 
netic calorimeter (ECAL) and anti-coincidence 
Veto counters. The superconducting magnet has 
the shape of a cylindrical shell with the inner 
diameter of 1.2 m and length of 0.8 m; it pro- 
vides a central dipole field of 0.8 Tesla. The eight 
layers of double-sided silicon tracker sensors are 
placed in planes transverse to the magnet axis. 
The silicon tracker measures the trajectory of rel- 
ativistic singly charged particles with an accu- 
racy of 10 fj, in the bending and 30 /i in the non- 
bending coordinates. It provides also measure- 
ments of the particle energy loss which allows to 
distinguish the charge. The timc-of-flight system 
(TOF) containing four detection layers, measures 
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singly charged particle transit times with an ac- 
curacy of 140 psec and also yields energy loss and 
coordinate measurements. The transition radi- 
ation detector (TRD) is situated on the top of 
the spectrometer and consists of twenty 12 mm 
thick foam radiator arrays, interleaved by arrays 
of 6 mm diameter gas proportional tubes filled 
with a Xe/C02 mixture. The TRD provides an 
e^/hadron separation better than one hundred 
up to an energy of 200 GeV as well as precise 
charged particle coordinate measurements. The 
RICH detector is installed below the last TOF 
plane and consists of a 3 cm thick aerogel radia- 
tor with a refraction index of 1.05, a mirror and 
pixel type matrix photo-tubes for the light detec- 
tion measures of the velocity of the single charged 
particle with an accuracy better than a fraction 
of a percent. The ECAL detector is situated at 
the bottom of the AMS02 setup. It is a three- 
dimensional (65x65xl7cm'^ electromagnetic sam- 
pling calorimeter with total length of 16Xq, con- 
sisting of 1 mm diameter scintillating fibers sand- 
wiched between grooved lead plates. 

3.2. Performance of photon detection 

Cosmic 7-rays may be detected in AMS by 
two different methods. The conversion mode in- 
volves the reconstruction in the tracker of the 
e'^e~ pairs produced by 7 conversions in the ma- 
terial upstream of the first layer of silicon sensors 
[33134135] . In the single photon mode, the 7-rays 
are detected in the electromagnetic calorimeter 

The performance of AMS02 detector for 7-rays 
has been studied with the AMS simulation and re- 
construction program based on GEANT j^Hl • The 
simulated performances have been validated us- 
ing the AMSOl data for the sub-detectors present 
during the shuttle test flight [321 and the test 
beam data obtained with prototypes of the new 
or modified modules. 

The Monte Carlo sample for the present study 
includes more than 10^ reconstructed events in- 
cluding both cosmic 7-rays and charged-particle 
backgrounds over the relevant energy range. The 
latter in order of decreasing importance include 
protons. He and C nuclei and electrons. With a 
charged particle background rejection of 0(10'*) 



to O(IO^), we obtain a background-to-signal ratio 
of the order of a few percent. The principal back- 
ground at this level is due to the galactic diffuse 
7-ray emission. 

3.2.1. The Conversion Mode 

The event signature for this mode are two re- 
constructed tracks in the Tracker originating from 
a vertex located upstream of the first silicon layer 
of the tracker'^. 

The incident 7-ray energy and direction were 
determined by adding the reconstructed momenta 
components of the pair, evaluated at the en- 
trance of the AMS02 detector. 

The main source of background are p and e~ 
which interact in the AMS detector, producing 
secondaries, mainly delta rays, which result in 
double-track events associated with a common 
origin at the interaction point. 

The conversion of the secondary photons pro- 
duced in the vicinity of the AMS, i.e. in the ISS 
body and solar panels, was found to be negligible 
in comparison to the expected 7-ray fluxes. 

The following criteria are applied to reject 
background events: 

• Identify events with interactions; 

• Identify charged particles entering the TRD 
from the top and fire all the tubes along its 
reconstructed trajectory. 

• Identify reconstructed large invariant mass 
events. 

• Identify particles entering the fiducial vol- 
ume of the AMS through the side of the 
TRD ). 

A preliminary rejection factor of 5 x 10^ was 
obtained for each different cosmic ray species (e~ 
and p) , after all selection cuts have been applied. 

3.2.2. The Single Photon Mode 

The event signature for this mode is the pres- 
ence of electromagnetic-type energy deposition in 

^The material in front of the first silicon tracker plane, 
consists of the TRD, the first two layers of ToF scintilla- 
tors, and mechanical supports, represents ~ 0.23Xo. 
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the ECAL, while almost nothing is found in the 
other AMS sub-detectors. 

The identified backgrounds contributing to the 
cosmic 7-ray signal are events with charged par- 
ticles^ either passing undetected in the gaps of 
the AMS active tracking volume or entering the 
ECAL from the side. 

The following criteria are applied to reject the 
background events: 

• Identify p, He by analyzing the 3- 
dimensional shower development in ECAL; 

• Identify charged particles by requiring the 
trajectory direction of the reconstructed 
ECAL shower passes inside the AMS sen- 
sitive volume and reject these events. 

The rejection factors for different cosmic ray 
species after all cuts have been applied are : > 
6x10" for e±, (2.5± 1) x 10^ for p and > 1.7 xlO^ 
for He nuclei. 

3.2.3. Acceptances and resolutions 

The simulations are used to parameterize the 
AMS performance for 7 detection in terms of ac- 
ceptance, effective area, angular and energy reso- 
lutions, and background rejection. Figure|21shows 
the acceptance and effective area for the two de- 
tection modes. The corresponding energy and an- 
gular resolutions are shown in Figure |3F- 

The parameterized performance is used to es- 
tablish the AMS02 sensitivity for the different 
scenarios in which high energy 7-rays are pro- 
duced by the annihilation of dark matter near 
the Galactic Centre. In a first approximation, we 
consider a 7-ray source located at galactic longi- 
tude /=0 and galactic latitude b=Q. 

3.3. The sensitivity to the gamma flux and 
Confidence Level determination 

We have developed a ROOT-based 40: sim- 
ulation program, the AMS-7 Fast Simulator 
(AMSFS) gH, in order to investigate the AMS 
capability to localize non-isotropic radiation, ei- 
ther point-like or diffuse. Here we describe the 

^mostly e~ , p and He nuclei 

^we have chosen a conservative estimate of the single pho- 
ton mode acceptance; a second study reports a 50% higher 
acceptance 1371 
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Figure 2. AMS acceptance as a function of 7-ray 
energy for the two detection modes (top panel). 
The effective areas versus zenith angle at 50 GeV 
7-ray energy (bottom panel). 



computational approach implemented in the sim- 
ulator. 

We define the detector's 7-ray source sensitiv- 
ity as the minimum flux required to achieve a 
specified level of detection significance. The sig- 
nificance S* of a detection is given by: 



S(E^ > Et) 



N°' 
7 



[E^ > Et) 



(7) 



where N°'"'{E^ 



^B{E^ > Et) 
> Et) and B{E^ > Et) are re 



■7 

spectively the total number of detected photons 
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from the source and the number of background 
photons faUing within the source area above an 
energy threshold Ef N°'"' and B are functions 
of the effective detection area A{E) of the instru- 
ment, the angular resolution expressed in terms 
of the solid angle i^{E), the observation time Tobs 
and the differential spectra: 



N;'%E^>Et)^ f 

J E 



, ^^A{E)Tobsdm)dE, 
Et Jn dE dil 

(8) 



and 



B{E^ >Et)= / ^-^A{E)T,bsdn{E)dE. 
JEt Jn dE dn 

(9) 
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Figure 3. The 68% containcment angular resolu- 
tion for both complementary detection modes as 
a function of energy (top panel). Energy resolu- 
tions as a function of original 7-ray energy (bot- 
tom panel). 



In order to establish the significance level of the 
observation, we require a minimum of three de- 
tected gamma events. 

We use the analytical expressions resulting 
from the best fit to the curves shown in Fig- 
ures 121 and for the energy dependence of the 
acceptance, the angular and energy resolutions. 
The solid angle over which the background is 
integrated for a given source is n{E)=TTalg{E), 
where ergs is the detector angular resolution de- 
fined within which 68% of the source photons fall. 

The calculation of the Galactic Center observa- 
tion time Toiis is based on a 3-year misssion on the 
International Space Station (ISS) g^.The AMS 
observation time is not uniformly distributed over 
the celestial sphere since the ISS is in a 51.6° or- 
bit, and the detector is fixed rigidly to the ISS. 
Taking into account the precession of the orbital 
plane of the station about the Earth's pole, a full 
sky coverage is obtained about 5.3 times per year. 
The exposure (effective area x the observation 
time) varies with the photon energy due to the 
energy dependence the effective area, and the po- 
sition in the sky, due to the orbit precession. 

The dependence of the effective area on the in- 
clination of the photon direction (0), the time 
dTobs spent by the detector viewing the Galactic 
Center within a specific viewing inclination d0, 
has been calculated and then integrated over the 
field-of-view (6* range up to 42° for the conver- 
sion mode and 22° for the single photon mode) 
and convoluted with the corresponding effective 
area: A{E,dd) x dTobs- The time intervals when 
ISS orbits over the South Atlantic Anomaly re- 
gion are excluded. 

The source spectrum dN^/dEdVL in equation © 
corresponds to the photon differential spectrum 
of the dark matter annihilation calculation incor- 
porating the halo profile model and the choice 
of the particle physics parameters, including the 
mass of the WIMP candidate. The background 
flux dB/dEdO. corresponds to the isotropic ex- 
tragalactic 7-ray background radiation and the 
galactic diffuse radiation (the latter is due mainly 
to the decay of tt^s produced by interactions of 
the cosmic rays with the interstellar medium). 
The extra-galactic component has been measured 
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by EGRET to be gS]: 



extragal. 

dE dVL 



where $o=(7.32 ± 0.34) x 10" 



(cm^ 



(10) 



s sr 



GeW)-\ ko=0.451 GeV and e=-2.10 ± 0.03. 

The galactic diffuse flux is enhanced toward the 
galactic center and the galactic disk as measured 
by EGRET. In our calculation we use the param- 
eterization of the differential flux provided in |22| : 



dB 



gal- 



dE dVL 



= rox 



E 



GeV)-\ (11) 



(cm^ s sr GcV) ^, ro=l 



where Fq^S.G x 
GeV and a^-2.7. 

Finally, the Galactic center point sensitivity 
has to take into account the profile of the energy 
spectrum of the photons produced as a function 
of the neutralino mass m^. For this purpose we 
define the following function: 



A(mx) 



Ig_^ME)dE 
L 



(12) 



A{m^) includes the weight of the detector accep- 
tance A{E) on the total rate of photons expected 
to be detected for a given differential flux 



A 

V 10 
10 
10 




from Tasitsiomi ei al. 



(GeV) 



Figure 4. The (N^av) as a function of with 
standard set of parameters as compared with pa- 
rameterization 44 . A threshold Et = 1 GeV has 
been assumed. A normalization factor of 2.4 was 
applied to obtain compatibility between the two 
calculations. This factor is used in our calcula- 
tions. 



dE 



1 dN-y (c 



(13) 



Therefore we can rewrite the function A as 



A(mx) 



'-i^AiE)dE 



/■OC 
Et 



OC dN^ 
~dE 



dE 



(14) 



where 



Et 



dN^ 
~dE 



dE = N^(Ey > Et) 



(15) 



is the total number of continuum 7-rays above en- 
ergy Et mainly due to the decay of tt*^ mesons pro- 
duced in jets from neutralino annihilation. For 



the computation of Nj we have considered the 
parametrization from ref. gl], after checking its 
compatibility with PYTHIA parameterization in- 
cluded in DSS as explained in the caption of Fig- 
ure |3| 



N^iE, >Et)^-l — 
6 Vm^ 



5 f Et 10 Et 



3 mv 



(16) 



The differential continuum spectrum, assuming 
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N^{Et = m^) = 0, is: 
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(17) 



With formulas (|16|l and (|17|l we can compute the 
function A(m^) for a given value of the neutralino 
mass Toy! 



(18) 



and obtain the corresponding Confidence Level. 

The total number of photons detected by AMS 
is defined as 



N°''iE^ > Et) - Tobs / —A{E)dE. (19) 

JEt "'^ 

or, using the A(m^) function 



Table 1 

The lightest neutralino mass m^, the niSUGRA 
parameters toq, tan/3, the relic neutralino densi- 
ties, i.e. ri^/i^ and the values {N^av) as described 
in the text in units of 10~^® cm^ s~^; (masses are 
in GcV and the stars indicate values from |47j'). 



model 


B 


G 


I 


K 


L 




98.3 


153.6 


143.0 


571.5 


187.2 


Too 


59 


116 


178 


999 


299 


tan(3 


10.0 


20.0 


35.0 


38.2 


47.0 




0.12 


0.12 


0.12 


0.11 


0.10 




0.12 


0.13 


0.13 


0.09 


0.10 


(n^av) 


1013 


1283 


8380 


29344 


33438 


(n^av)* 


782 


1032 


6303 


70903 


18739 



N;'"{E^> Et)^TobsHm^) -^dE. 



(20) 



According to the detectability criterion defined in 
eq. 0, the minimum detectable flux Fmin, which 
corresponds to a significance S{Ej > Et)=3, is 
derived by requiring Nj°''^ — 3 \/B. The defini- 
tion of: 



95 



d<^^ 



Et 



dE 



-dE = 



3 ^B{E^ > Et) 



(21) 



leads to conservative 95% — 99% Confidence Level 
values. 

4. Results 

4.1. mSUGRA and Benchmark point sim- 
ulations 

The DSS program provides values of the 7-ray 
fluxes for the SUSY benchmark models j45l46| 
and the so-called "wild scan" configurations of 
the mSUGRA parameters. 



The SUSY benchmark models have been pro- 
posed to provide a common way of comparing 
the SUSY discovery potential of the future ac- 
celerators such as LHC or Linear Colliders. The 
thirteen SUSY scenarios correspond to different 
configurations of the five mSUGRA parameters 
with the trilinear coupling parameter Aq set to 
0. The models fulfill the conditions imposed by 
LEP measurements, the — 2 result, and the 
relic density constraint of 0.094 < Q.^h'^ < 0.129. 
To derive the gamma-ray fiuxes for some of these 
benchmark models we use our current MC simu- 
lation programs: DSS, which was described pre- 
viously. In particular the value of fl^h^ is calcu- 
lated in the DarkSUSY part, while the simultane- 
ous use of the SUSPECT and DarkSUSY package 
allows to perform Renormalization Group Equa- 
tions evolution from the GUT scale to EWSB 
scale. 

Table presents values of the lightest neu- 



tralino mass 



^x ' 



the mSUGRA parameters toq. 



tort/3, the neutralino relic density Q~^h and the 
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Table 2 

The expected number of photons detected in 3 
years for different benchmark models and various 
dark matter halo profiles. Since the benchmark 
model flux values are low in the scanned ($-y, m^) 
plane, a 3 GeV energy cut threshold has been ap- 
plied for signal and diffuse gamma background 
calculations. The sensitivities {No- = number of 
standard deviations above diffuse gamma emis- 
sion) have been calculated only for photons de- 
tected in the tracker assuming 3.0 background 
photons. In the ECAL detector, we expect 196.0 
photons for the galactic diffuse gamma emission. 



model 


B 


G 


I 


K 


L 


AT NFW 

AT NFW 
1 cuspy 

AjAIoore 
7 


0.22 
9.2 
26.4 


0.14 
6.0. 
17.1 


0.94 
40.8 
117.3 


0.35 
15.2 
43.7 


2.48 
107.8 
309.9 


NFW 
cuspy 

AjAIoore 


0.04 
3.2 
8.9 


0.03 
1.9 
5.7 


0.18 
13.8 
39.6 


0.06 
4.7 
13.5 


0.45 
35.8 
102.8 



values {N^av) as described in section 3. 

The corresponding values of neutralino mass, 
tan(3 and mo are also quoted. A fine-tuning pro- 
cedure has been applied as in 0^1 (to fulfill the 
relic density constraints). The were obtained 
with fast simulation by the convolution of the dif- 
ferential 7-ray fiuxes with angular and energy res- 
olution, and applying the acceptance factors of 
the tracker (TR) and calorimeter (ECAL). 

The choice of the benchmark model sample was 
guided by the requirement of meaningful fiux val- 
ues. These values are low, however, in more fa- 
vorable astrophysical scenarios, the expected 
are enhanced by substantial factors varying from 
40 in case of the most cuspy NFW halo profile or 
about a hundred in case of a Moore profile. 



The results in Table Q are also compared to 
those in where a different mSUGRA Monte 
Carlo was used. Good agreement was found be- 
tween the results of the two calculations. The ob- 
served differences for (N^av), which are at most 
~ 25% between our results and those in may 
be explained by more refined interfacing. 

In Table El detected by AMS during 3-year 
observation and the significance values for the 
benchmark models are presented for the NFW 
halo profile with the standard set of parameter, 
the most cuspy NFW profile and the Moore pro- 
file. The diffuse 7-ray background has been eval- 
uated with the procedure described in section 3. 
The hadron contribution can be considered neg- 
ligible for a point-like source, as the proton sup- 
pression factors range between 10~^ and 10^^. 
For the standard NFW profile, only the values 
of the signal-to-background ratio are given as the 
expected values are not significant. 

4.2. Predictions from "wild scan" simula- 
tions 

4.2.1. mSUGRA results 

We have performed a "wild scan" in the 
mSUGRA parameter space. Six thousand 
models have been simulated in the region of 
0.0 < rJx^^ ^ 0-129. The values below the 
WMAP lower constraint on Q^h'^ (0.094) belong 
to the additional non-thermal neutralino produc- 
tion scenarios. 

The ranges of the mSUGRA parameters used in 
the simulation were: 

sign(/x) not constraint 
50. < mo < 3000. 
50. < mi/2 < 1600. 
0.1 < l^ol < 2000. 
3. < tan{P) < 60. 

The results for the integrated gamma fluxes from 
the Galactic Centre as a function of the Xi mass, 
presented in Figure [S] were obtained for a NFW- 
standard proflle, and for a 7-ray energy threshold 
oi Et = 1 GeV. Figure shows the results for the 
more favorable NFW cuspy dark matter proflle. 

4.2.2. AMSB results 

For the prediction of the gamma-ray flux in the 
AMSB framework, the scheme proposed by SUS- 
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PECT was used for the evolution of the AMSB 
parameters up to the EWSB, as for our mSUGRA 
simulations. Therefore we use the same DSS in- 
terface. 

The constraints set on the four AMSB param- 
eters (as described in 07|) were: 

sign(/i) not constraint 

10^ < M3/2 < 106. 

10^ < ao < 15 • 10^ 
3. < tan{l3) < 60. 

For completeness, we have checked that this ap- 
proach was compatible with the analytical ap- 
proximation as described below. 
From equation (2) we get: 



%{An,E^ > Et) = 1.19 10"^* TV. 



pb 



flTeV\' 



S19 cm s 



2„-l 



(22) 



and, under the hypothesis of a pure Wino LSP, 
the pair annihilation proceeds by exchange of a 
charged Wino. The neutral Wino, here assumed 
to be the WIMP-LSP, can annihilate into a W- 
boson pair (W^°1F° W'^W~). We have consid- 
ered the results of IHU^ for the parameterization 
of the corresponding annihilation cross section, in 
the non-relativistic limit: 



yrpb (23) 



where xw 



mi 



Thus the equation for the integral flux of pho- 
tons from Wino LSP annihilation becomes: 



4.3. Kaluza Klein results 

The main annihilation channels of LKP into 
Standard Model particles are charged lepton 
pairs for about 59% and quark pairs for about 
35% In our case, the calculation of sec- 

ondary gamma ray yield were based on the formu- 
las of cru(B(i)B(i)^//) [nj. The contribution 
of gamma rays produced from channels with lep- 
tons 48 has been neglected, thus providing more 
conservative results for the gamma fluxes at high 
energies. As in our previous calculations and fol- 
lowing we have used equation H22|l with again 
the value from equation H16|l . To obtain the 
flux for a given halo profile, we have used the cor- 
responding value of Sig given in the paragraph 



s 

o 

X 

-a 



NFW profile 




AMSB - DSS A AMSB - Analytic 
Kal Liza-Klein 

< n h- < 0.025 
0.025 < nh' < 0.094 

0.094 < fl h' < 0.129 (WMAP compatible) 
Benchmafk model points 



95-99 % C,L 



Neutralino or B^j^ Mass (GeV/c ) 



■ (2-xv^)2 ■ 



<l>^{An,E^ > Et) = 1.16 10"'^ iV. 

fTeVV ^ _2 _i 

X Ei9 cm s 

J 

(24) 

where we have used equation (|16|) for the iV^.^ value 
with the scaling factor as described in the caption 
of Figure^ 



Figure 5. The integrated 7 flux from the Galactic 
Centre as a function of for the NFW halo 
profile parameterizations with the standard set 
of parameters. The considered models are the 
mSUGRA scheme, AMSB scenario and Kaluza- 
Klein Universal Extra-dimensions. The various 
selections were done by varying flh^ cuts. 



-4 




Neutralino or B,,> Mass (GeV/c^) 



Figure 6. The integrated 7 flux from the Galac- 
tic Centre as a function of for a cuspy 
NFW dark matter halo profile as described in the 
text. The considered models are the mSUGRA 
scheme, AM SB scenario and Kaluza-Klein Uni- 
versal Extra-dimensions. The various selections 
were done by varying tan{[3) cuts. 



2.1, thus resulting in: 

^.t{m,E^ > Et) = 7.2 10"i^iV^x 

flTeVV 2 , (25) 

X Sig cm ^ 

\ m J 

We checked that our results are compatible 
with [H]. The expected 7 fluxes for the AMSB 
models and Kaluza-Klein models are also shown 
in Figures 5 and 6. 

4.4. Sensitivity for considered models 

The 95% CL was obtained by varying, within 
the uncertainties, the diffuse 7 background spec- 
trum as measured by EGRET 50^ in the Galactic 
Centre described in section 3.3. With 
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Figure 7. The resulting ^Li/H yield synthesized 
due to residual neutralino/Kaluza-Klcin annihi- 
lation during the epoch of Big Bang nucleosyn- 
thesis for the same models as considered before. 
Also shown is the central value of the ^Li/H as 
observed in the low-metallicity star HD84937. 



a required minimum of 3 photons over 3(T of 
the diffuse gamma background and a 1 GeV en- 
ergy threshold and 3 year exposure time, the 
sensitivity to Galactic Center measurements is 
(7.0 ± 0.4)10~^cm^^s^^ with only a small resid- 
ual dependence on m^. For both sets of astro- 
physical conditions, the predicted 7 fluxes for the 
AMSB and Kaluza-Klein models are above the 
95% CL for WIMP masses below 400 GeV. This 
indicates the potential of detection or exclusion of 
AMS02 in the case of the less conventional SUSY 
scenario with non-thermal production of neutrali- 
nos, or other dark matter candidate proposed by 
the Kaluza-Klein extra-dimension theories. 

It has been recently shown WD that residual 
dark matter annihilation during the epoch of Big 
Bang nucleosynthesis may result in an efficient 




U6 contribytion 
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production of ^Li. In Figure 7 we show the re- 
sulting ^Li/H ratio in the dark matter models 
studied in this paper. The ^Li/H yields have 
been calculated using the parametrisations given 
in Ref. jUl- 

The predicted abundances are compared to 
the value reported for the low-metallicity halo 
stars, such as HD84937, *^Li/H w 8.47 ± 3.10 x 
10~^^ (one of the first stars where a ^Li de- 
tection had been claimed). It is seen that, even 
in a context of a standard NFW profile (Fig- 
ure 5), the observed ®Li abundance is consistent 
with values produced with certain model config- 
urations, in particular in non-thermal scenarios 
for the AMSB model. Detections of ^Li have 
been reported for ~ 10 other stars |53I54| . with 
abundances comparable to IID84937. As this is 
far from what expected in cosmic ray scenarios 
which may synthesize ^Li, it is possible that the 
^Li abundance in low-metallicity stars is in fact an 
indirect signal of dark matter annihilation during 
Big Bang nucleosynthesis. 

Conclusions 

The DarkSUSY and SUSPECT programs were 
used to provide the 7-ray flux predictions from 
the Galactic Centre region, for the benchmark 
mSUGRA models, the AMSB scenario and 
Kaluza-Klein Universal Extra-dimensions mod- 
els, in order to evaluate the discovery potential of 
AMS for non-baryonic dark matter. Only mod- 
els of mSUGRA scenario with large tan(3 yield 
measurable signals on a realistic time scale. This 
conclusion is confirmed in a second study with a 
"wild scan" mSUGRA simulations. Various as- 
pects such as dependence of the results on heavy 
quark masses or CP-odd Higgs pole contributions 
may change the predictions |26j . The most sig- 
nificant signals correspond to the supersymmetric 
configurations with the neutralino mass ~ 100 
GeV. 

The sensitivity of the AMS-02 detector for 7- 
ray fluxes from a point-like source will allow 
to detect fluxes smaller by a factor of 2 to 3, 
compared to those measured by EGRET exper- 
iment in the GeV range, in the Galactic Cen- 
tre region [301 . AMS mission will also extend 



these measuremnets to a poorly explored energy 
range around 10 GeV, important for the detec- 
tion of a low mass neutralino. In the TeV range, 
the Galactic Centre was observed by ground- 
based Air Shower Cherenkov Telescopes ( ACTs) , 
which have detected several intensive astrophys- 
ical sources |55I56| . The signal from the cen- 
tral source, as observed by different ACT experi- 
ments, has been analyzed in the context of SUSY 
dark matter by [S7]. More recently, the HESS 
Collaboration has published a discovery of a dif- 
fuse 7 emission in the galactic plane in TeV range, 
nearby the SgrA* emitter, possibly produced by 
the hadronic interactions of the Galactic Cosmic 
Rays with a complex of molecular clouds j58j . 
These results indicate that the choice of the ob- 
served dark matter source is crucial, and the in- 
tense astrophysical environment may be an obsta- 
cle for an exotic signal detection. A more promis- 
ing type of dark matter source could be a nearby 
Dwarf Spheroidal Galaxy such as DRACO, Sagit- 
tarius or Canis Major, presenting lower standard 
astrophysical backgrounds. These sources have 
been considered by several authors for the SUSY 
dark matter flux predictions |59l60j , and their ob- 
servation campains have been scheduled by ACT 
telescope experiments in the coming future. 

In the frame of the present study with the 
Galactic Centre source, the 3 year observation 
with AMS detector would provide 95% CL ex- 
clusion limits for several mSUGRA models in 
the case of a favorable dark matter galactic halo 
configuration, such as the cuspy or very cuspy 
NFW profiles. Furthermore, if the halo is made 
of clumps with inner profiles of cuspy type, or 
if there is a strong accretion around the central 
black hole, the expected signal would increase 
by two orders of magnitude. For such cases, 
given the excellent energy resolution of the de- 
tector, the discovery of a dark matter annihi- 
lation signal would be possible. In particular, 
the non-thermal Susy Breaking scenarios, as in 
case of the AMSB model, result in cosmolog- 
ically significant ®Li abundances, which, when 
confronted with the results for ^Li abundances 
in low-metallicity stars offer interesting perspec- 
tives for indirect dark matter searches and the 
detection of an annihilation signal by AMS. We 
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conclude, that a survey of the Galactic Centre by 
AMS has the potential to contribute significantly 
to our understanding of dark matter. 

Acknowledgements 

We thank all our colleagues from AMS for ex- 
pressing their support and interest in this study. 

REFERENCES 



2. 
3. 



6. 



9. 

10. 

11. 



12. 



13. 



14. 



15. 



16. 



N.BahcaU and X.Fan, Astrophys. J. 504 
(1998) 1. 

N.Dalai and C.R.Keeton, |astro-ph/0312072 
M.Azzaro, F.Prada and C.M.Gutierrez, 

astro-ph/0310487, . 

D.Zaritsky et al. Astrophys. J. Suppl. 478 
(1997) 39. 

C.L.Bennett et al. Astrophys. J. Suppl. Ser. 
148 (2003) 1. 

G.Bertone, D. Hooper 
hep-ph/0404r75| Phys 



and J. Silk, 

Kept. 405 (2005) 

279. 

J.R. Primack, 2001, |astro-ph/0112255| To ap- 
pear in Proceedings of International School of 
Space Science 2001, ed. Aldo MorscUi (Fras- 
cati Physics Series). 

G. Jungman, M. Kamionkowski, K. Griest 
Phys.Rept. 267 (1996) 195-373. 
L.Bergstrom, Rep. Prog. Phys. 63 (2000) 793. 

H. Goldberg, Phys. Rev. Lett. 50 (1983) 1419; 
J. Ellis et al Phys. Lett. B127 (1983) 233. 
R. Barbieri, S Ferrara, C. A. Savoy, Phys. 
Lett. B 119 (1982) 343; A.H. Chamseddine, 
A. Arnowitt, P. Nath, Phys. Rev. Lett. 49 
(1982) 970; L.J. Hall, J. Lykken, S. Weinberg, 
Phys. Rev. D 27 (1983) 2359. N. Ohta, Prog. 
Theor. Phys. 70 (1983) 542. 

L.Randall, R.Sundrum, |hep-th/9810155| 
Nucl. Phys. B557 (1999) 79. 
G.Giudice et al., J. High Energy Phys. 9812 
(1998) 027. 

D. Hooper and L.Wang, Phys. Rev D69 (2004) 
035001. 

E. W. Kolb and R. Slansky, Phys. Lett. B 135 
(1984) 378. 

T.Appelquist, H.C.Cheng and B.A.Dobrescu, 
Phys. Rev. D 64 (2001) 035002. 



17. G.Servant, T.M.P.Tait, Nuclear Phys. B 650 
(2003) 391. 

18. G. Gentile et al., 2004MNRAS.351, 903. 

19. J.F.Navarro, C.S.Frenk and S.D.M.White, 
Astrophys.J. 462 (1996) 563-575. 

20. B.Moore et al, Phys.Rev. D64 (2001) 063508. 

21. P.Gondolo, J. Silk 1999, Phys. Rev. Lett. 83, 
1719. 

22. L.Bergstrom, P.UUio, J.H. Buckley, Astrop. 
Phys. 9 (1998) 137. 

23. P. UUio, H. S. Zhao, M. Kamionkowski, Phys. 
Rev. D64, (2001) 043504. 

24. J. Edsjo, et al., JCAP 0409 (2004) 004, 
astro-ph/0405414 

25. G. Bertone and D. Merritt, Phys. Rev. D72 
(2005) 103502. 

26. A.Falvard et al, 2002, astro- ph/0210184, As- 
tropart. Phys. 20 (2004) 467. 

27. H.P. NiUes Phys. Rep. 110 (1984) 1. 

28. A. Djouadi, J.-L. Kneur, G. Moultaka, 



hep-ph/0211331 



29. 'http: / / www. physt o . se/'edsjo /darksusy /' 

30. T.Moroi, L.RandaU, Nucl. Phys. B 570, 
(2000) 455-472. 

31. T.G.Rizzo, Phys. Rev. D 64 (2001) 
095010; C.Macesanu, C.D.McMuUen, 
S.Nandi, hep-ph/0201300, Phys. Rev. D66 

(2002) 015009; H.C.Cheng, K.T.Matchev, 
M.Schmaltz, Phys. Rev. D 66 (2002) 036005. 

32. The construction of the Alpha Magnetic 
Spectrometer (AMS) for the International 
Space Station, AMS Collaboration, to be 
published in Nucl. Instr. and Methods. 

33. R.Battiston et al., Astropart. Phys. 13, (2000) 
51. 

34. G.Lamanna, Nucl. Phys. B (Proc. Suppl.) 
113, (2002) 177. 

35. G.Lamanna, AMS-Note 2003-03-03 . 

36. V.Choutko, G.Lamanna, A.Malinin, Int. J. 
Mod. Phys. A, 17, (2002) 1817. 

37. L.Girard, Proceedings to Moriond - Astropar- 
ticle 2005 and references in. 

38. Geant4: A simulation Toolkit, Geant4 col- 
laboration, S. AgostineUi et al., NIM A 506 

(2003) 250. 

39. M. Aguilar et al., Phys. Rep.336 (2002) 331. 

40. http://root.cern.ch 

41. J.Bolmont, PhD Thesis, Univ. Montoellier II, 



16 



2005. 

42. I.Sevilla Noarbe , AMS-Note 2004-03-03. 

43. P. Sreekumaret al, Astroph. Journ. 494, 

(1998) 523. 

44. A.Tasitsiomi and A.V.Olinto, 
|astro-ph/0206040, Phys. Rev. D66 (2002) 
083006. 

45. M.Battaglia et al, Eur.Phys.J. C24 (2002) 
311-322. 

46. J.Ellis et al, Eur.Phys.J. C24 (2002) 311-322. 

47. P.Gondolo et al. , astro-pli/0406204 , JCAP 
0407 (2004) 008. ' " 

48. L. Bcrgstrom, et al, Phys. Rev. Lett. 94 
(2005) 131301. 

49.0 Bertone, G. Servant, G. Sigl , 

fhep-ph/0211342, Phys. Rev. D68 (2003) 
044008. 

50. S.D. Hunter et al., Astroph. Journ. 481, 
(1997) 205. 

51. K. Jedamzik, Phys. Rev. D 70 (2004) 063524; 
Phys. Rev. D 70 (2004) 083510. 

52. R. Cayrel et al, Astron & Astrophys. 343 

(1999) 923. 

53. M. Asplund, S. D'Odorico, P.E. Nissen, As- 
tron & Astrophys. 357 (2000) L49. 

54. P. E. Nissen et al, Astron & Astrophys. 348 
(1999) 211. 

55. F. Aharonian et al, HESS Collaboration, 
|astr o-ph/0408145 L. RoUand and J.Hinton, 
HESS collaboration, proceedings of ICRC 
conference, Pune, India (2005). 

56. M. Mori, CANGAROO collaboration, pro- 
ceedings of 28th ICRC conference (2003) vol 
5. 

57. S. Profumo, Phys. Rev. D72 (2005) 103521, 
astro-ph/0508628, 

58. F. Aharonian et al, HESS collaboration. Na- 
ture vol. 439 (2006) 439. 

59. N. W. Evans, F. Ferrer and S. Sarkar, 
astro-ph/0311145, 

60. D. Hooper et a/, astro-ph/ 03 1 1 1 50 



